Making Graphene Luminescent 
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We show that strong photoluminescence (PL) can be induced in single-layer graphene on using an 
oxygen plasma treatment. PL characteristics are spatially uniform across the flakes and connected 
to elastic scattering spectra distinctly different from those of gapless pristine graphene. Oxygen 
plasma can be used to selectively convert the topmost layer when multi-layer samples are treated. 



Graphene is at the center of a significant research 
effort [l]. Near-ballistic transport at room temperature 
and high mobility [i \A U] make it a potential material 
for nanoelectronics[5|, [y], especially for high frequency 
applications Q. Furthermore, its transparency and me- 
chanical properties are ideal for micro and nanome- 
chanical systems, thin-film transistors, transparent and 
conductive composites and electrodes and 
photonics [ll). There are two main avenues to modify 
the electronic structure of graphene. One is by cutting 
it into ribbons and quantum dots[l, @, [TH, [H, Q3], the 
other is by means of chemical or physical treatments with 
different gases .to reduce the connectivity of the tt elec- 
trons network flR fl6| . One of the most popular insu- 
lating chemical derivatives is graphene oxide (G0)(l6|. 
Bulk GO solutions and solids do also show a broad lu- 
minescence background P3, EH- Hydrogen plasma was 
used to controllably and reversibly modulate the elec- 
tronic properties of individual graphene flakes, turning 
them into insulators [15|. Aggressive ox yge n treatment 
was applied to create graphene islands [19|. However, 
thus far, no photoluminescence (PL) was seen from in- 
dividual graphene layers, either cut into ribbons or dots, 
or chemically treated, making graphene integration into 
optoelectronics still elusive. 

Graphene samples are produced by micro-cleavage of 
graphite on a silicon substrate covered with 100 nm 
Si02[l]. The number of layers is determined by a combi- 
nation of optical microscopy and Raman spectroscopy fioi 
l2lj |. Optical imaging at 473 and 514nm is done in an 
inverted confocal microscope. The beam is reflected 
by a splitter and focused with a high numerical aper- 
ture objective (NA=0.95). Raman spectra are measured 
at 514nm with a Renishaw micro- Raman spectrometer. 
The samples are then exposed to oxygen/argon (1:2) RF 
plasma (0.04mbar, 10W) for increasing time (1 to 6 sec- 
onds) . The structural and optical changes are monitored 
by Raman spectroscopy and elastic light scattering. PL 
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Figure 1: (a) Confocal PL image excited at 473nm (2.62eV) 
for a graphene sample oxidized for 3s. Scale bar 5/xm. The 
bright PL spots are spatially localized, (b) Uniform emission 
after 5s. Scale bar 10/xm. For position 3 in (b), PL is bleached 
intentionally by intense laser irradiation, (c) Spectra detected 
at the positions marked in (b). These have broad PL centered 
~700nm (1.77eV). (d) PL transients detected at the positions 
indicated in (b). The dynamics can be described by a triple- 
exponential with decay times ^40ps, 200ps and lOOOps. 



decay dynamics is recorded by time-correlated single pho- 
ton counting (TCSPC) upon pulsed excitation at 530 nm 
(2.34 eV), with a time-resolution of ±3ps. The acquisi- 
tion time per pixel is of the order of few tens ms. The 
spatial resolution is ~800nm. The power on the sample 
is well below lmW, to prevent photo-damage. 

Spatially resolved PL shows bright and localized emis- 
sion for short treatment times (Fig. [U(a)). For longer 
times, the PL is strong and spatially homogeneous 
(Fig.[T](b)), w ith a single broad band centered at ~700nm 
(L77eV), Fig. GJc). Intense laser excitation with power 
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Figure 2: Raman characterization as a function of treat- 
ment time, (a) Raman spectra; (b) D to G intensity ratio; (c) 
FWHM(G); (d) ratio of the slope of PL background (m) to 
the G peak intensity, 1(G) 



exceeding lmW leads to photo-bleaching and a PL max- 
imum blue-shift. The excited state decay dynamics of 
this photoluminecent graphene (PLG) is complex. The 
PL transients of Fig. [TJd) can be described by a three- 
exponential decay with lifetimes ~40, 200,1000ps, sub- 
stantially longer than those observed in semiconducting 
nanotubes and amorphous carbon 22, [23[. Remarkably, 
the PL transients are nearly uniform across the complete 
spectrum. This implies that spectral diffusion due to 
energy migration, typical for heterogeneously broadened 
systems, is absent (see Figure 6 in Methods). 

A Raman investigation gives further insights into the 
evolution from pristine graphene to PLG. Fig. 2 plots 
the Raman spectra and the main fitting parameters (see 
Methods). A broad PL background is evident in Fig 
2a for treatment times above 2s. This is quite different 
from the case of the hydrogen plasma treated samples 
of Ref . flH] . where no luminescence was observed. Fig 
2a also shows a significant increase of the D and D' in- 
tensities, and the D+D' combination mode ^2950 cm -1 , 
which requires a defect for its activation. Note that in 
defected graphene the relaxation of the backscattering 
condition results in significant broadening of the second 
order modes. Defect scattering also broadens the first or- 
der peaks, eventually merging G and D' in a single large 
G band for treatment times above Is. 

Fig 2b plots the evolution of the D to G peak in- 
tensity ratio, 1(D) /1(G). This first raises and then de- 
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Figure 3: White-light spectrum of PLG (black line) compared 
to pristine graphene (blue Line) [20]. The red line is a PLG 
model using a Cauchy function for the complex refractive in- 
dex (see text) (32l |. 



creases for increasing time. The D peak intensity is 
a measure of the number of defects [3, HH (see Meth- 
ods). Ref. [25|], noted that 1(D) /1(G) varies inversely 
with the cluster size L a in poly- and nano- crystalline 
graphites: I(D)/KG)=C(A)/L a , where C(514.5nm)~4.4 
nm from Refs.[25|, |26j]. This is known as the Tuin- 
stra and Koening relation (TK). TK holds until a criti- 
cal defect density. Since the D peak requires the pres- 
ence of sixfold rings, when the network starts losing 
them, 1(D) decreases with decreasing L a [24|]. In this case 
I(D)/I(G)=C'(A)L a 2 , with C'(514.5nm)-0.55nm- 2 [2i|. 
Combining the latter with TK, we deduce that Fig. 2b 
shows a continuous L a decrease down to^lnm, and a 
transition to a network with fewer sixfold rings for treat- 
ment longer than Is. This is further validated by con- 
sidering the evolution of the Full Width at Half Maxi- 
mum of the G peak, FWHM(G). In defect-free graphene, 
a variation of FWHM(G) is observed as a consequence 
of doping[13, 0, However, in the case of de- 

fected samples, peak broadening is a result of the activa- 
tion of q^O phonons. An empirical correlation between 
FWHM(G) and L a was reported in Ref.[30j considering 
a variety of disordered and amorphous carbons. Com- 
paring FWHM(G) in Fig.2c with the trend in Ref.Q, 
again we get L a ~lnm for the longest treatment. The 
large FWHM(G) also implies a distribution of L a around 
the average value. 

The ratio of the slope of the PL background (m), to 
1(G) is often used in disordered carbons as a measure of 
the PL strength, when comparing different samples|31]. 
We thus plot m/I(G) in Fig 2(d). This reaches a maxi- 
mum, then decreases for the l ong est treatment, consistent 
with the lack of PL in Ref.[19| after oxygen treatment 
targeted at layer removing. 

White-light scattering spectroscopy of PLG reveals 
clear differences to pristine graphene. Fig. [3] illustrates 
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Figure 4: Correlation between PL and layer thickness, (a) PL 
image; (b)scattering image of the same sample area.(c, 
d) Corresponding cross sections taken along the dashed lines 
in (a,b). PL is only observed from treated SLG, marked 1L 



the scattering spectra of both materials on SiC^/Si. 
While pristine graphene appears dark throughout the 
spectral range covered in the experiment, correspond- 
ing to a negative interferometric contrast [20|, PLG 
shows weaker contrast with a positive sign for wave- 
lengths smaller than 580nm. The spectrum of PLG 
yields the complex refractive index n' = A n + B n /\ 2 + 
i * (Ak + Bk/X 2 )j with Cauchy parameters A n =2.76, 
Afe=0.06, £> n =3000, £^ = 1500 for the graphene thickness 
of 0.34nm, comparable to those obtained for GO [H]. 

The data presented so far are taken for single 
layer graphene (SLG). A different behavior is observed 
for multi-layer graphene (MLG), which remains non- 
luminescent following treatment. Indeed, PL intensity 
and scattering contrast are directly correlated, as seen 
in Fig. [4] for flakes of different thickness. The scattering 
contrast for treated MLG does not scale linearly with the 
number of layers, as in the case of pristine SLG [20|. The 



MLG spectrum only features negative contrast, while the 
positive contribution below 580nm observed for PLG is 
absent. Scattering spectra from treated MLG can be rep- 
resented by a superposition of treated and pristine SLG. 

Oxygen plasma etching of graphite proceeds layer-by- 
layer [33j]. Thus, in our case only the topmost layer is 
affected. The absence of PL in MLG means that emis- 
sion from the topmost layer is quenched by subjacent 
un-treated layers. This opens the possibility of engineer- 
ing sandwiched hybrid structures consisting of PLG and 
a variable number of pristine graphene layers. 

Oxygen plasma etching is expected to yield CO and 
CO2, by successively removing carbon atoms. Etching of 
graphite occurs both in the basal plane and at defect s[34l|- 
The latter is consistent with our observation of point-like 
PL features for short treatment times (Fig. Ufa)). 

It would be tempting to interpret the PL emission as 
coming from electron confinement in sp 2 islands with 
an average size of ~lnm, as indicated by Raman spec- 
troscopy. Indeed, since electrons in graphene behave 
as massless particles, energy quantization due to con- 
finement is expected to open a gap SE w vph/2d w 
2eVnm/d. The resulting quantum confined energy for 
a quantum disk of diameter d=2nm is leV. The ob- 
served emission energy distribution translates into a di- 
ameter distribution ranging from 0.94 to 1.29nm, in 
agreement with the Raman estimation. In this case the 
large spectral width of the PL signal, ~0.5eV, could 
result from a superposition of overlapping bands with 
narrow linewidth centered at different size-controlled (or 
quantum confined) energies, corresponding to heteroge- 
neous broadening. Then, the optical properties of the 
PLG would resemble those of pi-conjugated polymer 
films, where a distribution of conjugation lengths trans- 
lates into a strong in-homogenously broadened density of 
states|35]. At room temperature, laser irradiation in the 
red would lead to selective excitation of a subset of quan- 
tum confined states. Then, spectral hole burning, i. e. 
the selective photobleaching of this subset of homoge- 
nously broadened lines, should be possible. This bleach- 
ing could be a photochemical modification or even a com- 
plete removal of the absorbing subset. As a result we 
would observe a spectral hole, i.e. the subset absorbing a 
certain color would not contribute to PL. Fig. 5 plots the 
ratio of PL measured before (Ipl) and after (I pl, bleached) 
exposure to high power (>600/iW) pulsed laser light at 
647nm. The PL is measured at 530nm for low excita- 
tion power (~ 10/iW). No spectral hole is observed in 
the detected spectral range, as would be expected for 
a heterogeneous ensemble of narrow bandwidth emitters. 
Instead, only an irreversible and uniform reduction of PL 
intensity occurs. For other bleaching energies in the red 
spectral range (760,800nm) the same uniform decrease is 
observed, while in the blue (473,514nm) the PL slightly 
shifts to shorter wavelength (see Fig.lc). 

Thus, we find that the observed large spectral width of 



4 




640 


i i 

680 


-I 1 L 

720 


760 


800 






X(nm) 






Kb) 1 




Figure 5: a) PL Intensity ratio for the area indicated by the 
circle in b). 



0.5eV reflects mainly homogenous broadening, uniform 
across the PLG sheet. This is supported by the ab- 
sence of spectral diffusion in the time-resolved data, ex- 
pected for heterogeneous films[35] (see Methods). If PL 
would indeed result from quantum confined states 
size-related heterogenoeus broadening would need to be 
far smaller, probably below O.leV, requiring a very nar- 
row size distribution of ~ ±0.04nm, instead of the ~ 
±0.18nm needed for the 0.5eV broadening. Since oxi- 
dation is expected to occur at different lattice sites and 
configurations, such high degree of ordering would seem 
unreasonable. Moreover, while for increasing oxidation 
times a successive decrease of the effective size distribu- 
tion would be expected, the spectral characteristics of 
the PL emission remain nearly constant. In conclusion, 
although the identification of L a as the quantum con- 
finement length of massless electrons would be tempting, 
we rather assign the observed PL to CO-related localized 
electronic states at the oxidation sites. 

In summary, we have shown that spatially uniform PL 
can be induced in single-layer graphene on substrates by 
selective plasma oxidiation. Remarkably, bi- and mutli- 
layer flakes remain non-luminescent, while their elastic 
scattering spectra indicate the formation of sandwich- 



like structures containing unetched layers. The resulting 
photoluminescent material could pave the way towards 
graphene-based optoelectronics. 
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Methods 
Raman Background 

Raman spectroscopy is a fast and non-destructive 
method for the characterization of carbons. These show 
common features in the 800-2000 cm -2 region: the G and 
D peaks. The G peak corresponds to the E23 phonon 
at the Brillouin zone centre. The D peak is due to the 
breathing modes of sp 2 rings and requires a defect for its 
It comes from TO phonqns around 

and is 



25 



activation [24 

K[i3, [111, is active by double resonance (DR)|36 
strongly dispersive with excitation energy due to a Kohn 
Anomaly at K[37|. The 2D peak is the second order of 
the D peak. This is a single band in monolayer graphene, 
whereas it splits in four in bi-layer graphene, reflecting 
the evolution of the band structure [21]. The 2D peak is 
always seen, even when no D peak is present, since no de- 
fects are required for the activation of two phonons with 
the same momentum, one backscattering from the other. 
DR can also happen as intra-valley process, i.e. connect- 
ing two points belonging to the same cone around K or 
K'. This gives rise to the D' peak,~1620cm -1 in defected 
graphite. The 2D' is the second order of the D' peak. 

Ref.[25|, noted that 1(D) /1(G) varies inversely with the 
cluster size L a in poly and nano crystalline graphites: 
I(D)/KG)=C(A)/L a , where C(514.5nm) - 4.4nm from 
Refs.[25l [5^1 . The original idea was to link 1(D) to 
phonon confinement. The intensity of the non- allowed 
D peak would be ruled by the defect-induced lifting of 
the Raman fundamental selection rule. Assuming that 
graphite becomes uniformly nano-crystalline, the D peak 
evolution can be estimated using Heisenberg indetermi- 
nation principle: 1(D) ex Aq, with AqAxocft and Ax~L a . 
We now know that the D peak activation is due to DR 
and not to phonon confinement. However, also in this 
case, the higher the number of defects, the higher the 
chance of phonon-defect scattering and, thus, the higher 
1(D). Again, since the G peak is not defect-activated, 
even within DR, one can expect TK to hold. Now L a 
is an average inter-defect distance, instead of a cluster 
size. This is a very simple picture, which has proven 
effective to compare graphitic samples for increasing dis- 
order. However, we note that a complete theory for the 
D and G Raman intensity and their dependence on the 
number of defects is still lacking. 
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Figure 6: PL transients measured at four different detection 
energies after excitation at 530 nm. These are nearly uni- 
form across the complete spectrum, indicating that spectral 
diffusion due to energy migration, typical for heterogeneously 
broadened systems, is absent. 



PL Transients 

Figure 6 plots the PL transients measured at four dif- 
ferent detection energies after excitation at 530 nm. All 
can be modeled by a tri-exponential decay with decay 
constants of 30, 250 and 1000 ps, with different relative 
contributions. Remarkably, there is no direct correla- 
tion between emission energy and decay dynamics, as 
would be expected for heterogeneously broadened sys- 
tems. In this case, spectral diffusion due to energy mi- 
gration would lead to faster decay in the blue spectral 
range and a delayed signal rise on the same time scale in 
the redjUj. The decay traces can also be modeled using 
a stretched-exponential model function. 
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